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Abstract: Fermion superpartners - neutralinos and charginos in the supersymmetric eco- 
nomical 3-3-1 model are studied. By imposition R parity, their masses and eigenstates 
are derived. Assuming that Bino-like is dark matter, its mass density is calculated. The 
cosmological dark matter density gives a bound on mass of LSP neutralino in the range of 
20 -7- 100 GeV, while the bound on mass of the lightest slepton is in the range of 60 -r- 130 
GeV 
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1. Introduction 

The Standard Model (SM) of high energy physics provides a remarkable successful descrip- 
tion of presently known phenomena. In spite of these successes, it fails to explain several 
fundamental issues like generation number puzzle, neutrino masses and oscillations, the 
origin of charge quantization, CP violation, etc. 

One of the simplest solutions to these problems is to enhance the SM symmetry 
SU(3)c SU(2)l U(l)y to SU(3)c «> SU(3)l U(l)x (called 3-3-1 for short) |, |, § 
gauge group. One of the main motivations to study this kind of models is an explanation in 
part of the generation number puzzle. In the 3-3-1 models, each generation is not anomaly 
free; and the model becomes anomaly free if one of quark families behaves differently from 
other two. Consequently, the number of generations is multiple of the color number. Com- 
bining with the QCD asymptotic freedom, the generation number has to be three. For the 
neutrino masses and oscillations, the electric charge quantization and CP violation issues 
in the 3-3-1 models, the interested readers can find in Refs. Q, |p and p, respectively. 

In one of the 3-3-1 models, the right-handed neutrinos are in bottom of the lepton 
triplets [^ and three Higgs triplets are required. It is worth noting that, there are two 
Higgs triplets with neutral components in the top and bottom. In the earlier version, these 
triplets can have vacuum expectation value (VEV) either on the top or in the bottom, but 
not in both. Assuming that all neutral components in the triplet can have VEVs, we are 
able to reduce number of triplets in the model to be two [0, Q (for a review, see P). Such 



a scalar sector is minimal, therefore it has been called the economical 3-3-1 model |10|. In 



a series of papers, we have developed and proved that this non-supersymmetric version is 
consistent, realistic and very rich in physics P, |l^, 11, 12]. 



In the other hands, due to the "no-go" theorem of Coleman-Mandula [^], the internal 
G and external P spacetime symmetries can only be trivially unified. In addition, the 
mere fact that the ratio Mp/Mw is so huge is already a powerful clue to the character 
of physics beyond the SM, because of the infamous hierarchy problem. In the framework 
of new symmetry called a supersymmetry [14, 15 1, the above mentioned problems can be 



solved. One of the intriguing features of super symmetric theories is that the Higgs spectrum 
(unfortunately, the only part of the SM is still not discovered) is quite constrained. 

It is known that the economical (non-supersymmetric) 3-3-1 model does not furnish 
any candidate for self-interaction dark matter [16| with the condition given by Spergel and 



Steinhardt |17|]. With a larger content of the scalar sector, the supersymmetric version 
is expected to have a candidate for the self-interaction dark matter. An supersymmetric 



version of the minimal version (without extra lepton) has been constructed in Ref. |18| 
and its scalar sector was studied in Ref. [|l9|. Lepton masses in framework of the above 
mentioned model was presented in Ref. [^], while potential discovery of supersymmetric 
particles was studied in [^. In Ref. |22], the R - parity violating interaction was applied 
for instability of the proton. 

The supersymmetric version of the 3-3-1 model with right-handed neutrinos 0] has 



already been constructed in Ref. j^Sj. The scalar sector was considered in Ref. |24| and 



neutrino mass was studied in Ref. |25]. Note that there is three- family versions in which 



lepton families are treated differently pM and their supersymmetric versions are presented 



in Ref. |27]. It is worth mentioning that in the previous papers on supersymmetric version 
of the 3-3-1 models, the main attention was given to the gauge boson, lepton mass and Higgs 
sectors. An supersymmetric version of the economical 3-3-1 model has been constructed 
in Ref. |^^. Some interesting features such as Higgs bosons with masses equal to that of 
the gauge bosons - the W and the bileptons X and Y, have been pointed out in Ref. p9|. 
Sfermions have been considered in Ref. [pQ ]. 

In a supersymmetric extension of the (beyond) SM, each of the known fundamental 
particles must be in either a chiral or gauge supermultiplet and have a superpartner with 
spin differing by 1/2 unit. Both gauge and scalar bosons have spin-2 superpartners with 
the electric charges equal to that of their originals: called neutralinos without electric 
charge and charginos if carrying the latter one. In the Minimal Supersymmetric Standard 
Model (MSSM), in some scenario, the neutralino can be the lightest and plays a role of 
dark matter. In this paper, we will focus an attention to neutralinos and charginos in the 
supersymmetric economical 3-3-1 model. 

This article is organized as follows. In Sec. we present fermion and scalar content 
in the supersymmetric economical 3-3-1 model. The necessary parts of Lagrangian is also 
given. In Section y, we deal with neutralinos sector. To find eigenstates and their masses, 
we have to adopt some assumptions. Section |^ is devoted for charginos. In Section ^ we 
present analysis of relic neutralino dark matter mass density and the limit on its mass. 
Finally, we summarize our results and make conclusions in the last section - Sec. |6|. 



2. A review of the model 

In this section we first recapitulate the basic elements of the supersymmetric economical 
3-3-1 model |28|. R — parity and some constraints on the couplings are also presented. 



2.1 Particle content 

The superfield content in this paper is defined in a standard way as follows 

F = iF,F), S = {S,S), V = iX,V), (2.1) 

where the components F, S and V stand for the fermion, scalar and vector fields while 
their superpartners are denoted as F, S and A, respectively |14, p^]. 



The superfield content in the considering model with an anomaly-free fermionic content 
transforms under the 3-3-1 gauge group as 

QiL = [ui, di, u'j ~ (3,3,1/3), 
uIl, S'£~(3M,-2/3), d5i~(3*,l,l/3), 

QaL = \da,-Ua,d'a] ~(3,3*,0), a = 2, 3, 

uIl ~ (3*, 1, -2/3) , d^L, €l ~ (3*, 1, 1/3) , 

where the values in the parentheses denote quantum numbers based on (SU(3)c, 811(3)1,, 
U(l)x) symmetry. i?£ = {vrY and a = 1, 2, 3 is a generation index. The primes superscript 
on usual quark types (n' with the electric charge q^i = 2/3 and d' with q^' = —1/3) indicate 
that those quarks are exotic ones. 

The two superfields x and p are at least introduced to span the scalar sector of the 
economical 3-3-1 model [10|: 



X= (x?,r,X^)'^~ (1,3,-1/3), 



p=(p+,pO,p+)^~ (1,3,2/3). 

To cancel the chiral anomalies of higgsino sector, the two extra superfields x' ^-iid /? must 
be added as follows 

x'=(x'i°,x'+,x'2°)^~(l,3M/3), 
p'=(p'f,p'0,p'2")^~(l,3*,-2/3). 

In this model, the SU(3)l (X" U(1)x gauge group is broken via two steps: 

SU(3)l U(l)x ^ SU(2)i ^ U(l)y '''"-^"' U(1)q, (2.2) 



where the VEVs are defined by 






(0,r;,0), V2{pY ={0,v',0). 



(2.3) 



The VEVs w and w' are responsible for the first step of the symmetry breaking while u, u' 
and V, v' are for the second one. Therefore, they have to satisfy the constraints: 



u, u' , V, v' <C w, w' . 



(2.4) 



It is emphasized that the VEV structure in ( p.3| ) is not only the key to reduce Higgs 
sector but also the reason for complicated mixing among gauge, Higgs bosons, etc. As it 
will be shown in the following, the mentioned VEV structure causes flavour violation in 
the D-term contributions. 

The vector superfields Vc, V and V containing the usual gauge bosons are, respectively, 
associated with the SU(3)c SU(3)l and U(l)x group factors. The colour and flavour vector 
superfields have expansions in the Gell-Mann matrix bases T" = A"/2 (a = 1,2, ...,8) as 
follows 



Vr. 



1 



A"y„ 



Vr 



1 



A"*Fe„; V = h'^Va, V = -lx''*Va 



where an overbar indicates complex conjugation. For the vector superfield associated 
with U(l)x, we normalize as follows 



XV' = {XT^)B, T^ 



1 

76 



diag(l,l,l). 



The gluons are denoted by 5° and their respective gluino partners by A^, with a = 1, ... ,8. 
In the electroweak sector, V"" and B stand for the SU(3)l and U(l)x gauge bosons with 
their gaugino partners Xy and A^, respectively. 

With the superfields as given, the full Lagrangian is defined by Cgusy + ^soft, where 
the first term is super symmetric part, whereas the last term breaks explicitly the super- 
symmetry [^] . The interested reader can find more details on this Lagrangian in the above 
mentioned article. In the following, only terms relevant to our calculations are displayed. 

2.2 i?-parity 

For the further analysis, it is convenience to introduce i?-parity in the model. Following 
Ref. 1 25], i?-parity can be expressed as follows 



R - parity = (-l)25(_i)3(B+£) 
where invariant charges C and B (for details, see Ref. [^) are given by [^ 



(2.5) 
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3. The neutralinos sector 

The higginos and electroweak gauginos mix each with other due to effects of the electroweak 
symmetry breaking. The neutral higginos and gauginos combine to make the mass eigen- 
vectors called neutralinos. In this section, the mass spectrum and mixing of the neutralinos 
is considered. 

The gauginos mass terms come directly from the soft term given by 



Csoft = Y, MbW'w' + M^BB. 



(3.1) 



6=1 



Because of the R-parity conservation, the higginos mixing terms come from the //—term 
determined as 



^/ I ^^/ 



-fi-term = PxXX + fJ-pPP ■ 



(3.2) 



Finally, the mixing terms between higginos and gauginos are a result of Higgs-higginos- 
gauginos couplings 



C = -V2g ((/>*T» A'^ - V2g\+^ (V'+T^ 



(3.3) 



Expanding Eqs ( |3.lD , (3.2) and ( |3.3| ), we obtain the neutralino mass matrix in the gauge- 

eigenatates basis ■0° = ixiiXi 1X2^X2 ^Pi^ Pi t^^^^^^&t'^^''^*]^ which is given in the 
Lagrangian form 



c= [r] M^r 



(3.4) 



with the following notations 



~ ^ W4 + iWs ~ ^ W4 - iW5 



(3.5) 
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(3.6) 



where M^ = M^ = M45. The mass matrix M^ can be diagonahzed by an unitary matrix 
U to get the mass eigenstates. It means that we can find matrix U satisfying: 



UMU 



-1 



"Ni ' ''^N2 ' "^A^s "^Af4 ' ""^Afs ' "^A^e ' ^Nr ' 



(3.7) 



Diag(m^^ 

with real positive entries on the diagonal. 

In general, the parameters Mb,M3,Ms, M45, fi^, fip can take arbitrary complex phase. 
However we can choose a convention to make MB,M3,M8,Mi5 to be all real and posi- 
tive. If we choose the parameter f^^, fJ-p to be real and positive then we must pick up the 
ix) ) ix') ) (p) 1 {p') to be real and positive too. If ^j^. and /ip are not real, then we obtain 
the CP violating effects in the potential. Therefore, as the same as in the MSSM [15|, it 
is convinience to choose the Mx' A*p to be real but without fixing the sign of ^j^, /ip. 

Getting exact eigenvalues and eigenstates of the mixing mass matrix (p]B) is very 
difficult task. Hence, we make some assumptions which is suitable for theoretical comments; 
and their correctness could be tested by the future experiments. 

In this paper, we assume that 



v,v' ,u,u' ,w,w' <C l/ip 



Mb\ , l/ip - A^3| , l/^p - Ms\ ,\lJip-M 



45 1 



and 



v,v' ,u,u' ,w,w' < \^y, - Mb\ ^^lx --^sl ,1% --^sl ,IMx --^45! 



(3.^ 



(3.9) 



In the above limit, using a small perturbation on the neutralinos mass matrix (p.6|), we can 
obtain the neutralino mass eigenstates, which are nearly a "higginos-like" , a "Bino-like", 
a "zino-like", an "extrazino-like ", a "xino-like", and the conjugated of the "xino-like" 
corresponding to 



7Vi=i3,7V2 = >V3,iV3 = W8,iV4 

' - p"±pT ~ 



X\N. = X, 



Nq,N7 



V2 



.N8,Nq 



jVio,ivii = 



V2 



V2 



(3.10) 



with the mass eigenvalues: 



m 



'Ni 



Mb + 



g 



.12 



{u + u'f + {w + w'f g'^ {u-u'f + {w-w'f 



^N. = -^3 + 



108 {Mb + f^x) 

g'Hv-vf g'Hv + vf 
27 {Mb - fip) 27 {Mb + fip) 
g^ {u + u' 



+ 



108 {Mb - l^x) 



^ 27 {Mb -i^p)^ 



gHv + vf gHv-v'f 

i{M3-i^p) 8 {Ms + Up)' 



^N, = -^8 + 



{u-u'f + A{w-w'f 
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24(7W8 + %) 
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m 



24(7W8 + Mp) 24(A^8-/"p)' 

45 



_ g2 ^2n^uu' + 7W45 (^^^ + u'^)] 



m 



o2 ^2n^yj^' + 7W45 (u;2 + w'^)] 



where 



"^iVe = 



"^N, = 



"^TVs = 



'^N,= 



"^iVio = 






m-all = 



I gHv-v'f gHv-v'f g'Hv + v'f 

^'^ 8{f,p- M3) 24 {fip - Ms) 27 {fip - Mb) 

I g^v + vf gHv + vf g'Hv-vf 
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/^xl + 2 "^«ll + "^"22 - Y (w-all - "la22)^ + 4m2^2 

/^xl + ^ 



"Ifell + "lfe22 - Y ("Ifcll - "1622)^ + 4m2^2 



"lall + "1-022 + Y ("T-all " "la22)^ + 4777.2^2 



/^xl + 



"Ifell + "1-622 + Y ("Ifell - "1-622)^ + 4m^^2 



126 



-^^— ^ - + 9g^ 
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Mb-^Ix V%--^3 /Ix--^ 



(u + u'"] 



./\2 
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(3.11) 
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(3.12) 



We emphasize that A^g, A^3, A^8)-^45 were taken real and positive and ^^x-ii^p ^^^ ^^^^ 
with arbitrary sign. The mass values depend on the numerical values of the parameters. 
In particular case, we assume M.b < -^3 < -^8 < -A^45 *^ f^x'l^p- ^^ *^^^ case, we obtain 
the neutralino lightest supersymmetric particle (LSP), which is a Bino-like A'^i. In the 
following, we will focus our attention to the neutralino LSP. 

4. The charginos sector 

The charged winos (W"'", W~,3^^,3^^) mix with the charged higginos {x^ , x'^i Pi^, P2^ , 
Pi~, P2~) to form the eigenstates with the electric charges ±1. They are called charginos. 
As the same as in the MSSM, we will denote the charginos eigenstates by Q . The entries 
of the elements in the charginos mass matrix come from ( |3.l| ) , (|3.2| ) and ( ^.31) . In the gauge- 
eigenstate basis ip^ = (W+, y'^ , pi^, p2^ , ^^ ■, W^, y~ ^ P\~ ■, P2~ , X~), the chargino 
mass terms in the Lagrangian form are given by 



C, 



charginomass 



i;^] M~i;^ + H.c 



(4.1) 



with the Mr having the 2x2 block form: 
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^^~=U- 



(4.2) 



where A4 is 5 x 5 matrix given by 
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(4.3) 



In principle, the mixing matrix for positive charged left-handed fermions and negative 
charged left-handed fermions are different. Therefore, we can find two unitary 5x5 matrices 
U and V to relate the gauge eigenstates with the mass eigenstates 



(4.4) 



This means that the charginos mass matrix can be diagonalized by two unitary matrices 
U and V to obtain mass eigenvalues 
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(4.5) 
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To finish this section, we note that in the mo 


del under consideration there are five charginos; 


and they are subject of the future sti 


rdies. 









5. Neutralino dark matter 

In the model under consideration, there are eleven neutralinos Nn (n = 1, 
them is a linear combination of eleven R = —1 Majorana fermions, i.e. 



., 11), each of 



N„ 






(5.1) 



where A^^^ are the normalized eigenvectors of the neutralino mass matrix (3.6). The question 
to be addressed is that our consideration below comes with the conditions ( |3.8| ), ( |3.9D and 
Mb < -M-s < M.8 < -M-Ab ^ P'x^l^p- Assuming that the neutralino LSP is a Bino-like 
A^i, we should show its predicted relic density is consistent with the observational data. 
To answer the question, we must calculate cross section for neutralino annihilation and 
compare it with the observational data on dark matter by the WMAP experiment [p2[ 



nnMh' = (0.12771°:^^) 



(0.02229 ±0.00073). 



(5.2) 



In (^), the normalized Hubble expansion rate h = 0.73_qq3. We adopt the allowed region 
as 

0.0895 < nDAih"^ < 0.1214. (5.3) 

Before calculating, we should note that a precise determinations of the relic density 
requires the solution of the Boltzmann equation governing the evolution of the number 
density n 



"iV 



dn 
IE 



-3—n — iva) in 
a 



n 



eql 



(5.4) 



with a is the cross section of the A'j's annihilation and v is the relative velocity. The 
thermal average {va) is defined in the usual manner as any other thermodynamic quantity. 
In the early Universe, the species Ni were initially in thermal equilibrium, n^ = ricj^q. 
When their typical interaction rate T^ became less than Hubble parameter, F^ < H, the 
annihilation process froze out. Sine then their number in comoving volume has remained 
basically constant 



For the present purpose, we will use approximate solution for Xf = -^ 



x/=ln 



my I 45 i 



(5.5) 



where 5* stands for the effective energy degrees of freedom at the freeze-out temperature 
(y/^ ~ 9) and Gn is the Newton constant. Typically one finds that the freeze-out point 
Xf is basically very small (~ 20)' "^^^ relic mass density p^ at the present is given in |33| 



with the suppression factor ( -tS- I ~ ^ following from the entropy conservation in a 
comoving volume. The coefficients a and b are determined by 



/ 



"Af 



1 p 



^9{m^- ruf) --^ [{A} + BJ) {4m% - m)) + QAfBfm)\ (5.7) 



where p = W (M- — m?r] and At and Br will be defined below. The sum is taken over 



■ I't) anu j-Vf ana d t 

the different types of particle-antiparticle pairs into which the A^ annihilate. 

In order to calculate the LSP mass density, to determine the Aj and Bj coefficients, we 
need to write down the low-energy effective Lagrangian from interactions. The calculation 
of the annihilation cross section in our model is straightforward in principle but quite 
complicate in practice. To ease our work, we consider only the most important channels 
for neutralino annihilation in the lowest order (tree-level) of perturbation theory for the case 
in which the LSP is a nearly pure Bino A'^i. The most important channels are annihilation 
into a pair of fermions 

NiNi^ff,{f = q,l,u) (5.8) 

and into a pair of charged Higgs scalar 

NiNi^H+H~,H^H^. (5.9) 

Because the Bino does not couple to W^ , Z and Z\ there is no annihilation of pure 
Bino to W+W- and ZZ, Z' Z or to Z' Z' . 



10 



Now we list the couplings needed in computation of the annihilation cross sections. 
The couplings of Bino B to quarks and leptons and their two scalar partners are given by 
the following piece of Lagrangian: 



13_ 
V3 



- (llb - llb) + (rrB - frB 



1 






1 



^QiQi - ^u^ui + ^d^di - ^u^u'" + ^dpd'l ) B 

The couplings of neutral Higgs and charged Higgs are determined in the following 
terms 



^/3 
2 
3 



-3 [xxB - XXBJ + - [XX'B - XX B 



+ '^{ppB--p~pB]-^[~pp'B--p'~p'B 



2 ^~i j; 
3 



(5.10) 



With the help of the mentioned couplings, the Feynman diagrams for Bino annihilation 
processes are depicted in Fig. || 
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Figure 1: Feynman diagrams contributing to annihilation of Bino dark matter 



We note that the LSP can annihilate to the particles only if theirs mass is lighter than 
the LSP mass. In the [^], by studying the Higgs sector, we have obtained one charged 
Higgs with mass equal to the W-gauge bosons mass {mw) and the other ones have mass 
equal to the bilepton mass My > 440 GeV. Therefore, in the region mf% < mw, the LSP 
cannot annihilate to charged Higgs and the top-quark as well as the exotic quarks and only 
the annihilation channels into ordinary fermion pairs such as A^A^ — > //, except for the 
top-quark, are available. 

From the Feynman diagram for Bino annihilation processes, the effective Lagrangian 
for a Majorana fermion N interacts with an ordinary quark or lepton / can be written 
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down: 

^e// = J] ^7^75^77/. (AfPL + BjPr) f (5.11) 

/ 



with 



YfJ- YfJ- 



12m1 12m1 ' 
/-L fn 

^2 J2 v2 J2 



Yf g'^ Yf g'- 

Bf = ^4 -^ (5-12) 

^ 12m?- 12m?- ^ ' 

where 1/,, Y/j are hypercharge of left- and right-handed ordinary quark and lepton. 

In deahng with Eq.(^.6|), we have taken into account g' = 0.6 in the model under con- 
sideration and suggested that all squarks mass are heavier than all sleptons and especially, 
m-q = 5mj. In Fig. B the LSP mass density dependence on its mass has been plotted 



From Eqs. (p^), (|5.7]) and (5.12), it follows that the density increases for increasing 



of sfermion mass (oc m?) and decreasing of the LSP mass (oc ^^). Fig. ^ shows also that 
the LSP mass is in the range of 100 GeV. 

In Fig. ^, the LSP mass density dependence on two dimensional space of parameters 
LSP mass and sparticle mass has been plotted. The LSP density is drawn as plane. We 
have divided the space of parameters into allowed and disallowed regions, where boundaries 
of acceptable region according to (^]^) are drawn as grid green plane and grid blue plane. 
From the Fig. ^ we obtain the lighter sfermion mass is heavier than Bino mass. We also 
obtain the bounds for mass of the sfermions: 60 GeV < mj < 130 GeV, while the masses 
of the LSP is in the range of: 20 GeV < m^ < 100 GeV. It should be noted that this 
result coincides with estimation given in ||3^ (see Fig 1 in page 1114). 

Let us consider the case m^ = mj. The LSP mass density has been plotted in Fig. 
(H). The figure shows that the LSP mass density is very small; it is even smaller than the 



lower bound given by the |32]. This means that this case is excluded by the WMAP data. 

6. Conclusions 

In this paper we have investigated the neutralinos and charginos sector in the supersym- 
metric economical 3-3-1 model. Accepting conversational assumption such as in the MSSM, 
eigenmasses and eigenstates in the neutralinos sector were derived. By some circumstance, 
the LSP is Bino-like state. 

In the charginos sector, the mass matrix can be diagonalized by two 5x5 matrices V 
and U . 

Assuming that Bino-like is dark matter, its mass density is calculated. 

The cosmological dark matter density gives a bound on mass of LSP neutralino is in 
the range of 20 -^ 100 GeV. In addition we have also got a bound on sfermion masses to 
be: 60 -^ 130 GeV. We have also shown that the case ms = mr is excluded by the recent 
experimental WMAP data. Our result is favored the present bound and it should be more 
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Figure 2: 

LSP's mass density as a function of its mass. The blue, red, yellow, green, violet curves 

are allowed by 

and lower experimental limits given in 



nij = 50,60,100,160 GeV, respectively . The horizontal lines are upper 



cleared in the near future. As in the MSSM, the neutralinos in our model gain the masses 
in the working region of the LHC. Consequently they could be checked in coming years. 
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